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Abstract 

Epitaxial strain is a proven route to enhancing the properties of complex oxides, however, 
the details of how the atomic structure accommodates strain are poorly understood due to the 
difficulty of measuring the oxygen positions in thin films. We present a general methodology for 
determining the atomic structure of strained oxide films via x-ray diffraction, which we demonstrate 
using LaNiOa films. The oxygen octahedral rotations and distortions have been quantified by 
comparing the intensities of half-order Bragg peaks, arising from the two unit cell periodicity of the 
octahedral rotations, with the calculated structure factor. Combining ab initio density functional 
calculations with these experimental results, we determine systematically how strain modifies the 
atomic structure of this functional oxide. 
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Strain engineering is an appealing route to controlling electronic [I], ferroic [2H1], and 
structural properties in complex oxide heterostructures [5]. The sensitivity to strain 
arises from the strong electron-, spin-, and orbital-lattice coupling inherent to transition 
metal oxides [7H9]. In ABO3 perovskites (Fig. [T^,), the flexible corner-shared -BC-6 octahedral 
network allows epitaxial strain to be accommodated by either bond-length distortions (Fig. 
[IJ},c), and/or by changes in the magnitude and phase of the octahedral rotations (Fig. [Tji,e). 
Previous work has demonstrated how strain-based manipulation of bond-length distortions can 
be used to dramatically enhance ferroelectricity [TU1HT]. A less explored, yet perhaps equally 
productive [T51 H3] strategy is the control of bond angles to manipulate magnetism, metal- 
insulator transitions, and superconductivity through tuning of superexchange pathways and 
electronic bandwidths. However, the current understanding of strain-bond angle relationships 
is limited, in part due to the experimental challenge of measuring oxygen positions in thin 
films [TH - flB] . This limited understanding of how strain alters both bond lengths and angles 
is a clear obstacle to engineering tailored functionalities into oxide films and heterostructures. 

In perovskite oxides, octahedral rotations offset the oxygen atoms from the face-centered 
positions effectively doubling the pseudocubic unit cell, producing a distinctive set of half- 
order Bragg peaks depending on the tilt pattern [T7j. The octahedra can rotate about the 
pseudocubic [1 0] (a), [0 1 0] (6), and/or [0 1] (c) directions [T8H20] ; these rotation angles 
are referred to as a, (3, and 7, respectively, and are depicted in Fig. [T^,. To describe the 
phase of the octahedral rotations along each axis, we employ the Glazer notation, in which 
a superscript is appended to each axis to indicate whether neighboring octahedra rotate 
in-phase ( + ), out-of-phase (~) or if rotations are absent (°) |18J . In the absence of cation-site 
displacements, the allowed and forbidden half-order diffraction peaks are a direct fingerprint 
of the octahedral rotations in perovskite oxides. 

While half-order peaks have been employed previously to identify crystallographic space 
groups of perovskite films [21], in this work we demonstrate how these Bragg peaks can be 
used to quantify the octahedral rotations and bond angles of strained perovskite oxide films. 
The oxygen positions, and thus the bond angles and lengths, are obtained by comparing 
the intensities of a systematic series of half-order Bragg peaks with the calculated structure 
factor of the oxygen octahedra. Combining the diffraction results with density functional 
calculations, we examine the evolution of the rotation angles with strain for a model functional 
oxide. The relative simplicity of the experimental approach, requiring only high flux x-ray 



2 




Figure 1: The available structural distortions afforded in perovskite oxides through heteroepi- 
taxial thin film growth, a, Rotations of the oxygen octahedra are pervasive in bulk perovskites 
and commonly decomposed into components a, /3 and 7 about each of the Cartesian axes. In 
coherently strained perovskite films, the BOq octahedra can distort through a contraction (b) 
or elongation (c) of the in-plane B-0 bond lengths (d) under compressive or tensile strain, 
respectively. Simultaneously, the octahedra can accommodate the substrate-induced in-plane 
lattice parameters by rotating perpendicular to the substrate (d), and/or rotating along the 
axes parallel to the substrate plane (e). 

diffraction and the calculation of structure factors, makes this a promising strategy for 
determining the full atomic structure of a broad range of strained perovskite films. 

We chose to study LaNiC>3 (LNO) as an example because it exhibits octahedral rotations 
along each Cartesian axis. Bulk LNO has an a~a~a~ rotation pattern (a = /3 = 7 = 5.2°) 
with a R3c symmetry and a pseudocubic lattice parameter of 3.838 A [22]. This results 
in uniform Ni-O-Ni bond angles of 165.2° and Ni-0 bond lengths of 1.935 A. The isotropic 
nature of bulk LNO is ideally suited to exploring how strain alters bond lengths and angles in 
an anisotropic manner. Additionally, LNO has garnered interest due to its novel properties 
such as electric- field control of conductivity [23], correlated Fermi liquid behavior [21], and 
recent predictions of superconductivity [25] . 

Our epitaxial LNO films were deposited on SrTiOs (a = 3.905 A) and LaA103 (a = 
3.795 A) substrates by molecular beam epitaxy [22] in order to probe how both tensile (1.7 
% on STO) and compressive (-1.1% on LAO) strain are accommodated by changes in the 
bond angles and lengths in the perovskite structure. The films grown on SrTi03 (STO) and 
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LaAlC>3 (LAO) are 173 A and 95 A thick, respectively, and coherently strained with c-axis 
parameters of 3.807 and 3.895 A as confirmed through reciprocal space maps and (00L) scans. 
All x-ray diffraction measurements were performed at room temperature. 

In addition to exploring how strain modifies the octahedral distortions, growing on the 
two substrates, STO and LAO, enables us to investigate how the substrates' structural 
properties interact with the octahedral rotations of the LNO films through hetero epitaxy. 
For example, STO is cubic at room temperature and lacks octahedral rotations, whereas 
LAO is rhombohedral with the a~a~a~ rotation pattern. Thus, by measuring the structural 
distortions in the LNO/STO films, we can determine whether octahedral rotations are 
suppressed when growing on a cubic substrate. Similar analysis of the distortions in the 
LNO films on LAO enables us to determine if the rotational pattern present in the substrate 
couples to that found in the film. 

The rotation pattern of the LNO films is identified by observing the presence and absence 
of specific half-order Bragg peaks. We find half-order peaks in the films grown on both 
substrates (Fig. [2^,,b), which confirms that heteroepitaxial growth on a cubic substrate 
does not suppress octahedral rotations by forcing atomic registry between the film and 
the substrate; the epitaxial relationship is not simply "cube-on-cube." Both films exhibit 
half-order peaks when the following conditions are met: 1) h, k, and I are equal to n/2, where 
n is an odd integer and 2) h — k ^ I, k — I ^ h, and h = I ^ k. Peaks are not observed 
where h = k = I, such as ( 1 /2 x /2 1/2), or where one of the indices is equal to an integer, such 
as ( 1 /2 1 3 /2). This reciprocal lattice indicates the presence of an a~a~c~ tilt pattern within 
both films [T7j, with a = (5 due to the biaxial strain state of the films (a = b). 

The intensities of the half-order peaks can be compared to calculated peak intensities 
in order to quantify the octahedral rotation angles. The calculated peak intensities are 
determined from Eq. Q: 

1 1 4 

I = h ■ f-s ■ (0 ^ {Y,D 3 \F hkl \ 2 ), (1) 
sm(?7) sm(29) ^—^ 

where Iq is the incident photon flux, l/sin(?7) is a correction for the beam footprint, l/sin(2@) 
is the Lorentz polarization correction, Fhki is the structure factor for the oxygen atoms, 
and Dj is the relative volume fraction of structural domains discussed below. The incident 
photon angle is 77, while the scattering angle is O. The structure factor for each Bragg peak 



4 





(1/2 1/2 3/2 




- (3/2 1/2 1/2) 




■jlJ 





0.4 0.5 0.6 1.5 1.6 
L (STO r.l.u) 



10 6 




|lao Ilao 




LNO 




10 5 




S i 








10 4 




f (3/2 3/2 5/2/ 1 


10 3 




A 



1.4 1.5 2.4 2.5 
L (LAO r.l.u) 

Peak Mens. I/I Calc. I/I 



(11!) 


1.0 


1.0 


[222) 


1.2 


1.0 


hi* 

(222) 


0.46 


0.49 


119 

{2 2 2) 


0.04 


0.08 


(222) 


5.8 


6.0 


(2 2 2) 


1.3 


1.3 



Peak 


Meas. I/I 


Calc. I/I„ 


fi T si 


1.0 


1.0 


[ill] 


2.4 


2.2 


(HI) 


0.56 


0.60 


[111) 


1.2 


1.3 


(11!) 


0.09 


0.20 


(!!!) 


44.1 


50.5 


(!!!) 


39.9 


31.7 


(!!!) 


18.2 


20.0 


(!!!) 


37.9 


25.7 


(222) 

(!!!) 


3.0 
0.75 


3.6 
0.51 


(!!!) 


0.62 


0.61 



Figure 2: Half-order peaks of LNO films grown on STO and LAO. a, In LNO on STO, the 
h = k diffraction peaks are more intense than those with h ^ k, as a > 7. b, The opposite 
intensity pattern is observed for the LNO on LAO, as a < 7. Intense half-order peaks 
from octahedral rotations within the LAO substrate are also present; this intensity from the 
substrate is less than 10% of the film peak intensity, c, Comparison of the measured and 
calculated peak intensities for LNO on STO for a = (3 = 7.1° and 7 = 0.3°. d, Comparison 
of the peak intensities for LNO on LAO for a — (3 — 1.7° and 7 = 7.9°, with a domain ratio 
of 1:2.15:0:0. 

is obtained from, 

24 

Fhki = fo 2 - 2J exp[27ri(/iw n + kv n + lw n )], (2) 

n=l 

where the position of the n th oxygen atom within the unit cell is given by (u, v, w), and 
fo 2 - is the O 2- form factor [27] . The La and Ni atoms are assumed to reside on the ideal 
corner and body-centered positions, respectively, and thus we do not consider a rhombohedral 
distortion of the La positions. Given this assumption, only the oxygen atoms contribute 
to the half-order peak intensities. We account for the octahedral rotations by doubling the 
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unit cell in each direction, giving rise to a non-primitive unit cell consisting of 24 oxygen 
positions (8 octahedra). The oxygen positions are displaced from the ideal face-centered 
positions by the appropriate rotation angle. To calculate the intensity of a given peak, the 
relative volume of all possible structural domains (Dj) is needed. When calculating the 
structure factor, these domains determine whether the octahedron nearest to the chosen 
origin rotates clockwise or counterclockwise about each axis. Within this tilt pattern there 
are four domains that give rise to different structure factors: one where a, (3, and 7 are 
all positive, and three where two of the angles are positive and the third is negative. The 
domain volume fractions are obtained from a set of symmetrically equivalent half-order peaks 
with fixed I, for instance (rtr 1 ^ ± 1 /2 2/2), the relative intensities of which depend directly on 
the domain populations. We employ standard non-linear regression analysis to determine 
the optimal values of a and 7 and their respective errors. 

The films grown on cubic STO and rhombohedral LAO substrates exhibit different domain 
populations. In the LNO/STO film, we find an equal volume fraction of the four possible 
domains, indicating that the four domains are energetically equivalent when LNO is grown 
on a cubic substrate. In contrast the LNO/LAO film shows unequal domain populations 
between only two domains with an approximate 1:2 ratio. To explain this difference, we note 
that our LAO substrates contain extensive twinning, which we observed as two distinct peaks 
in cu-scans rocked about a (0 /) peak. For our LAO substrate the ratio of the two peaks 
is also 1:2 suggesting that the octahedral rotations and domain pattern of the substrate 
are transferred to the film. Since the rhombohedral distortion determines the trigonal axis 
about which the octahedra rotate in bulk LAO, this result demonstrates that substrates 
may be used to texture octahedral rotations by choosing appropriate substrate orientations. 
Heteroepitaxy therefore affords an additional mechanism by which epitaxial thin film growth 
can be used to engineer the physical properties of complex oxides. Additionally, imprinting 
of rotational behavior may play an important role in the interfacial properties of oxides. 

For LNO on STO, we obtain a = (3 = 7.1 ± 0.2° and 7 = 0.3 ± 0.7°. The measured and 
calculated intensities are given in Fig. pp. For LNO on LAO, we obtain a = (3 = 1.7 ± 0.2° 
and 7 = 7.9 ± 0.9°; the peak intensities are provided in Fig. |2jl. Knowing the tilt angles, and 
thus the oxygen positions, the atomic structures are obtained (Fig. [3]). 

As anticipated, the epitaxial strain strongly modifies the in-plane (equatorial) Ni-0 bond 
lengths (d^i-ci) due to coherency across the heterointerface. The in-plane djsa-o decreases 
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systematically from LNO/STO (1.968 ± 0.002 A) to bulk LNO (1.935 A) to LNO/LAO 
(1.916 ± 0.005 A) as the a- and 6-axis parameters of the substrate are reduced. The measured 
out-of-plane (apical) d^i-o-, however, are less sensitive to the strain; we obtain 1.933 ± 0.002 
A and 1.950 ± 0.002 A for LNO/STO and LNO/LAO, respectively. (Supplementary Table 
SI provides a full summary of our experimentally determined bond lengths and angles.) 
Interestingly we find that bi-axial strain only weakly modifies the in-plane Ni-O-Ni bond 
angle (0) of 165.8 ± 0.5° and 164.0 ± 2.0° for LNO/STO and LNO/LAO, respectively. In 
contrast, along the out-of-plane direction 9 is highly sensitive to the substrate- induced strain; 
this bond angle decreases from 175.2 ± 0.6° in LNO/LAO to 159.9 ± 0.6° in LNO/STO. We 
note that both the a (or 0) and 7 angles contribute to the in-plane bond angles and since 
we find that empirically the average of these angles is robust to strain, the insensitivity of 
the equatorial bond angle is not surprising, albeit counterintuitive. 

To investigate the origin of these effects and to determine the change in octahedral rotations 
for intermediate strain values, we perform density functional theory (DFT) calculations within 
the local-spin density approximation (LSDA) plus Hubbard U method ({7 e fr = U—J = 3 eV) as 
implemented in the Vienna ab initio Simulation Package (vasp) |28fl3~T] . We first optimize the 
internal degrees of freedom for bulk rhombohedral LNO (10-atom unit cell) at the experimental 
lattice parameters to confirm that we accurately describe the electronic and atomic structure. 
Consistent with the experimental data, our optimized structure is indeed metallic and consists 
of equal and alternating phase rotations of the NiOg octahedra about each Cartesian axis 
with a = (3 = 7 = 5.76°, in good agreement with the 1.5 K structural measurement of 5.37° 
[22] . We next carry out structural optimizations for the two homoepitaxially strained films 
using our experimentally determined lattice parameters. The substrates are not explicitly 
simulated in our calculations, but the symmetry reduction imposed by enforcing equal 
in-plane lattice parameters on the bulk R3c structure is included by optimizing the internal 
degrees of freedom within a monoclinic (20-atom) unit cell (space group C2/c). The excellent 
consistency between our experimental and theoretical structures is apparent in the structural 
data obtained for the nickelate films grown on STO and LAO as shown in Fig. [3^. To quantify 
the comparison between the two methods, we use a symmetry adapted mode analysis |32| 133] 
and take scalar products between the DFT and experimentally determined vectors £ that 
contain the full irreducible representations of the atomic displacement patterns of the strained 
nickelate films about the bulk rhombohedral structure. We find that £(dft) • £(exp) ~ 0.99 
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Figure 3: Atomic structure of LNO grown on STO and LAO as determined by x-ray diffraction 
(Experiment) and density functional theory (Theory). In a), the top data set corresponds to 
LNO/STO, while the bottom data set corresponds to LNO/LAO. The space group is C2/c 
(#15), with a = 9.406 A, b = c = 5.523 A, and a monoclinic angle of 125.95° for LNO/STO 
and a = 9.460 A, b = c = 5.367 A, and a monoclinic angle of 124.57° for LNO/LAO. The 
relation of the monoclinic unit cell (a,b,c) to the ideal perovskite (a p ,b p ,c p ) is a p = c/a/2, b p = 
b/y/2, and c p = y/6a. b, Tensile strain from the STO substrate reduces the 7 rotation, while 
increasing the a and rotations, c, Compressive strain induced from the LAO substrate has 
the opposite effect: increasing the 7 rotation, while decreasing the a and (5 rotations. 



for LNO on both LAO and STO which indicates excellent accuracy between the two structure 
determination methods. We note that the DFT ground state indicates an additional energy 
lowering distortion through small antiparallel displacements of ~ 0.01 A for the La atoms 
along the [110] directions. Such displacements would produce weak but measurable ( x /2 l /2 1 /2) 
and ( 3 /2 3 /2 3 /2) peaks, which however were not observed in the diffraction measurements, and 
likely persists in the DFT calculations due to the common exchange- correlation functional 
underestimation of the equilibrium volume. 
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We next extend our analysis to intermediate strain values to investigate the transition 
region between the experimentally obtained strain values. In this section, we carry out 
structural relaxations about the theoretical LSDA+f/ equilibrium volume rather than at 
the experimental lattice parameters reported previously. We explore bi-axial strain states 
by fixing the in-plane lattice constants at each strain value, and performing full structural 
optimizations of the internal atomic coordinates and c-axis parameter about the LSDA+[7 
reference structure. In Fig. [4^, we show that over the strain range investigated the NiOg 
octahedra always rotate out-of-phase about each axis with the a~a~c~ rotation pattern 
with only changes in their magnitude about each axis. In agreement with the experimental 
data, we find that compressive (tensile) strain predominately induces rotations of the NiC>6 
octahedra about the c-axis (in-plane axes). We attribute this effect to the epitaxial constraint 
that forces coherency between the film and the substrate, which contracts (elongates) the 
equatorial bond lengths and elongates (contracts) the apical bonds under compressive (tensile) 
strain (Fig. |3Jd). The a~ 'a~ c~ octahedral rotation pattern is able to accommodate the bond 
length distortions by increasing the NiOg rotation about the c-axis at the expense of the 
in-plane rotations. For the tensile case the amount of in-plane rotation decreases accordingly 
to sustain the bond length elongation in-plane. We find that there is an additional structural 
phase that is approximately 2 meV lower in energy for small tensile strains and shows an 
intriguing small charge disproportionation (CDP). The deviations in rotation angles between 
the two phases shown in Fig. [4^, are negligible, and since our room temperature experiments 
are likely unable to access this phase, we do not discuss it further. 

The transition region which captures the crossover from mainly rotations about the c-axis 
to the in-plane axes occurs near 0.50% strain or when the axial ratio c/a approaches unity 
(Fig. [lb). At this strain value, standard crystallographic distortion parameters that measure 
deviations in ideal six-fold coordination through the apical and equatorial bond lengths and 
the Ni-O-Ni bonds all show a sharp discontinuity and indicate almost ideal octahedra |34J. 
In the absence of the small CDP phase, we identify this abrupt, first-order phase transition 
to be is o symmetric [35j since the atomic structure remains monoclinic — the full symmetry of 
the C2/c space group is retained and the same Wyckoff positions are occupied — across the 
transition. The origin for first-order phase transition about c/a = 1 can be understood on 
geometric grounds (Supplementary Fig. S5): when the ideal NiOg octahedra is recovered at 
0.5% strain, the bi-axial lattice distortion imposed on the film by the substrate leads to a 
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Figure 4: The DFT calculated rotation angles a and 7 (a), Ni-0 bond lengths (b) and 
the c/a axial ratio (c) for the monoclinic nickelate structure as a function of strain (open 
symbols). A sharp first-order transition induced by heteroepitaxy occurs at approximately 
c/a = 1. The solid symbols (broken lines) in (a) indicate the bond angles associated with a 
lower symmetry structure described in the text. 

tiling of corner-connected octahedra that is incompatible with size of the unit cell. Recently a 
strain-induced isosymmetric phase transition, with the axial ratio being the order parameter 
has also been reported in insulating thin films of the rhombohedral antiferromagnetic and 
ferroelectric BiFeC>3 j6]. Since LNO and BiFeOs have distinctly different physical properties 
(LNO is metallic and non-polar), we suggest this intriguing strain-induced isosymmetric 
phase transition may be a universal feature of epitaxially strained rhombohedral oxide films 
with the a" a~ c~ rotation pattern. In principle such strain-induced changes to the rotation 
patterns can be used to tailor the bandwidth anisotropy, and thus macroscopic properties 
such as conductivity and magnetism in the thin film. In the case of the large bandwidth 
compound LNO, however, our DFT calculations (Supplementary Fig. S6) and experiments 
indicate that the strain induced by the STO and LAO substrates is insufficient to produce 
dramatic changes in the electronic structure. 

In summary, we have demonstrated a general strategy for quantifying oxygen octahedral 
rotations in strained complex oxide films. We have determined the atomic structure of 
strained LaNiOs, in which the NiOe octahedra are distorted and rotated to accommodate 
strain. It is anticipated that the future application of this method will enhance fundamental 
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studies of strain accommodation and relaxation in perovskite oxides, allow for more direct 
comparisons between bulk and thin film materials, and help elucidate the physical mechanisms 
of novel properties at oxide interfaces. 

Methods 

X-ray diffraction measurements were performed at the Advanced Photon Source on Sectors 
33-BM and 6-ID using four-circle and six-circle Huber diffractometers. When analyzing the 
peak intensities, we choose to use the peak height as a measure of intensity. Scans along 
h and I confirm that the peak widths are independent of their location on the reciprocal 
lattice. The measured intensity is determined by normalizing for the incident intensity and 
subtracting background intensity. The calculated intensity is determined using Eqns. (1) and 
(2). The corrections for the beam footprint, Lorentz polarization factor, and the q dependence 
of the O 2- form factor are taken into account in the calculated intensity. The adjustable 
parameters are the rotation angles, a and 7. In the case of LNO/LAO, we did not attempt 
to remove any diffuse background from the LAO substrate, which exhibits strong half-order 
peaks. We estimate the magnitude of this substrate-induced background to be less than 
10% of the intensity of the film peak. This value is less than the average agreement between 
the measured and calculated peak intensities. As such, we do not believe this introduces a 
significant error into our determination of a and 7. 

We carry out first-principles density functional theory (DFT) calculations within the 
local-spin density approximation (LSDA) plus Hubbard U method [28] as implemented in the 
Vienna ab initio Simulation Package (VASP) [30l [31] . We follow the Dudarev approach [29] 
and include an effective Hubbard term U e f / = U — J of 3 eV to accurately treat the correlated 
Ni 3c? orbitals. The core and valence electrons are treated with the projector augmented 
wave (PAW) method |36j with the following valence electron configurations: 5p 6 5rf 1 6s 2 
(La), 3p 6 3c? 9 4s 1 (Ni), and 2s 2 2p 4 (O). The Brillouin zone integrations are performed with 
a Gaussian smearing of 0.05 eV over a 5x5x5 Monkhorst-Pack fc-point mesh [37] centered 
at r, and a 450 eV plane-wave cutoff. For structural relaxations we also used a Gaussian 
broadening technique of 0.05 eV and relaxed the ions until the Hellmann-Feynman forces are 
less than 1 meV A -1 . In all calculations ferromagnetic spin order is enforced, and for our 
choice of t/ e //, we accurately reproduce the enhanced paramagnetic spin moment |38j . 



11 



Acknowledgements 

Work supported by UChicago Argonne, LLC, Operator of Argonne National Laboratory 
(Argonne). Work at Argonne and use the Advanced Photon Source at Argonne was supported 
by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, 
under Contract No. DE-AC02-06CH11357. We acknowledge use of the oxide-MBE at the 
Center for Nanoscale Materials. We acknowledge support from NDSEG (JMR) the NSF 
under grant no. DMR 0940420 (NAS) and discussions with the Army Research Office MURI 
program (award no. W911-NF-09-1-0398). Portions of this work made use of the SGI Altix 
COBALT system at the National Center for Supercomputing Applications under grant no. 
TG-DMR-050002S, and the CNSI Computer Facilities at UC Santa Barbara under NSF 
award no. CHE-0321368. Correspondence regarding the diffraction measurements should be 
made to P.R. (pryan@aps.anl.gov). Correspondence regarding the DFT calculations should 
be made to J.M.R. (rondo@mrl.ucsb.edu). 



* Electronic address: smay@ drexel.edu| 
[1] Ward, T. Z., Budai, J. D., Gai, Z., Tischler, J. Z., Yin, L. & Shen, J. Elastically driven 
anisotropic percolation in electronic phase-separated manganites. Nature Phys. 5, 885-888 
(2009). 

[2] Tokura, Y. & Nagaosa, N. Orbital physics in transition-metal oxides. Science 288, 462-468 
(2000). 

[3] Schlom, D. G., Chen, L. Q., Eom, C. B., Rabe, K. M., Streiffer, S. K. & Triscone, J. M. Strain 
tuning of ferroelectric thin films. Ann. Rev. Mater. Res. 37, 589-626 (2007). 

[4] Ramesh, R. & Spaldin, N. A. Multiferroics: Progress and prospects in thin films. Nature Mater. 
6, 21-29 (2007). 

[5] He, F., Wells, B. O., Ban, Z.-G, Alpay, S.-R, Grenier, S. M., Shapiro, S. M., Weidong, S., 
Clark, A. & Xi, X. X. Structural phase transitions in epitaxial perovskite films. Phys. Rev. B 
70, 235405 (2004). 

[6] Zeches, R. J., Rossell, M. D., Zhang, J. X., Hatt, A. J., He, Q., Yang, C.-H., Kumar, A., Wang, 
C. H., Melville, A., Adamo, C, et al.. A strain-driven morphotropic phase boundary in BiFe03. 



12 



Science 326, 977-980 (2009). 
[7] Torrance, J. B., Lacorre, P., Nazzal, A. I., AnsaldoE. J. & Niedermayer, C. Systematic 

study of insulator-metal transitions in perovskites i?Ni0 3 (i?=Pr,Nd,Sm,Eu) due to closing of 

charge-transfer gap. Phys. Rev. B 45, 8209-8212 (1992). 
[8] Chmaissem, O., Dabrowski, B., Kolesnik, S., Mais, J., Brown, D. E., Kruk, R., Prior, P., Pyles, 

B. & Jorgensen, J. D. Relationship between structural parameters and the Neel temperature in 

Sri^Ca^MnOs (0 < x < 1) and Sri^BayMnOs (y < 0.2). Phys. Rev. B 64, 134412 (2001). 
[9] Mizokawa, T., Khomskii, D. I. & Sawatzky, G. A. Interplay between orbital ordering and 

lattice distortions in LaMn0 3 , YV0 3 , and YTi0 3 . Phys. Rev. B 60, 7309-7313 (1999). 
[10] Choi, K. J., Biegalski, M., Li, Y. L., Sharan, A., Schubert, J., Uecker, R., Reiche, P., Chen, 

Y. B., Pan, X. Q., Gopalan, V., Chen, L.-Q., Schlom, D. G. &: Eom, C. B. Enhancement of 

ferroelectricity in strained BaTi0 3 thin films. Science 306, 1005-1009 (2004). 
[11] Haeni, J. H., Irvin, P., Chang, W., Uecker, R., Reiche, P., Li, Y. L., Choudhury, S., Tian, 

W., Hawley, M. E., Craigo, B., et al., Room-temperature ferroelectricity in strained SrTi0 3 . 

Nature 430, 758-761 (2004). 
[12] Zayak, A. T., Huang, X., Neaton, J. B. &: Rabe, K. M. Structural, electronic, and magnetic 

properties of SrRu0 3 under epitaxial strain. Phys. Rev. B 74, 094104 (2006). 
[13] Hatt, A. J. k, Spaldin, N. A. Competition between polar distortions and octahedral rotations 

in epitaxially strained LaA10 3 . (2008), URL |arXiv : 0808 . 3792 (unpublished) 
[14] Xie, C. K., Budnick, J. I., Hines, W. A., Wells, B. O. & Woicik, J. C, Strain-induced change 

in local structure and its effect on the ferromagnetic properties of Lao.sSro.5Co0 3 thin films. 

Appl. Phys. Lett. 93, 182507 (2008). 
[15] Woodward, D. I. &: Reaney, I. M., Electron diffraction of tilted perovskites. Acta Cryst. B 61, 

387-399 (2005). 

[16] Jia, C. L., Mi, S. B., Faley, M., Poppe, U., Schubert, J. & Urban, K., Oxygen octahedron 
reconstruction in the SrTi0 3 /LaA10 3 heterointerfaces investigated using aberration-corrected 
ultrahigh-resolution transmission electron microscopy. Phys. Rev. B 79, 081405(R) (2009). 

[17] Glazer, A. M. Simple ways of determining perovskite structures. Acta Cryst. A 31, 756-762 
(1975). 

[18] Glazer, A. M. The classification of tilted octahedra in perovskites. Acta Cryst. B 28, 3384-3392 
(1972). 

13 



[19] Woodward, P. M. Octahedral tilting in perovskites. I. Geometrical considerations. Acta Cryst. 
B 53, 32-43 (1997). 

[20] Woodward, P. M. Octahedral tilting in perovskites. II. Structure stabilizing forces. Acta Cryst. 
B 53, 44-66 (1997). 

[21] He, F., Wells, B. O. & Shapiro, S. M. Strain phase diagram and domain orientation in SrTi03 

thin films. Phys. Rev. Lett. 94, 176101 (2005). 
[22] Garcfa-Muhoz, J. L., Rodriguez-Carvajal, J., Lacorre, P. & Torrance, J. B. Neutron-diffraction 

study of i?Ni03 (i?=La,Pr,Nd,Sm): Electronically induced structural changes across the 

metal-insulator transition. Phys. Rev. B 46, 4414-4425 (1992). 
[23] Scherwitzl, R., Zubko, P., Lichtensteiger, C. & Triscone, J.-M. Electric-field tuning of the 

metal-insulator transition in ultrathin films of LaNi03. Appl. Phys. Lett. 95, 222114 (2009). 
[24] Eguchi, R., Chainani, A., Taguchi, M., Matsunami, M., Ishida, Y., Horiba, K., Senba, Y., 

Ohashi, H. & Shin, S. Fermi surfaces, electron-hole asymmetry, and correlation kink in a 

three-dimensional Fermi liquid LaNi0 3 . Phys. Rev. B 79, 115122 (2009). 
[25] Hansmann, P., Yang, X., Toschi, A., Khaliullin, G., Andersen, O. K. & Held, K. Turning a 

nickelate Fermi surface into a cupratelike one through heterostructuring. Phys. Rev. Lett. 103, 

016401 (2009). 

[26] May, S. J., Santos, T. S. & Bhattacharya, A. Onset of metallic behavior in strained 
(LaNi0 3 )„/(SrMn0 3 )2 superlattices. Phys. Rev. B 79, 115127 (2009). 

[27] Hovestreydt, E. On the atomic scattering factor for 2 ~. Acta Cryst. A 39, 268-269 (1983). 

[28] Anisimov, V. I., Aryasetiawan, F. & Liechtenstein, A. I. First-principles calculations of the 
electronic structure and spectra of strongly correlated systems: the LDA+U method. J. Phys.: 
Condens. Matter 9, 767-808 (1997). 

[29] Dudarev, S. L., Botton, G. A., Savrasov, S. Y., Humphreys, C. J. & Sutton, A. P. Electron- 
energy-loss spectra and the structural stability of nickel oxide: An LSDA+U study. Phys. Rev. 
B 57, 1505-1509 (1998). 

[30] Kresse, G. & Furthmiiller, J. Efficient iterative schemes for ab initio total-energy calculations 

using a plane-wave basis set. Phys. Rev. B 54, 11169-11186 (1996). 
[31] Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave 

method. Phys. Rev. B 59, 1758-1775 (1999). 
[32] Stokes, H. & Hatch, D. Isotropy Subgroups of the 230 Crystallographic Space Groups (World 

14 



Scientific, Singapore, 1988). 
[33] Orobengoa, D., Capillas, C, Aroyo, M. I. & Perez-Mato, J. M. AMPLIMODES: Symmetry 

mode analysis on the Bilbao crystallographic server. J. Appl. Cryst. 42, 820-833 (2009). 
[34] Robinson, K., Gibbs, G. V. & Ribbe, P. H. Quadratic elongation: A quantitative measure of 

distortion in coordination polyhedra. Science 172, 567-570 (1971). 
[35] Christy, A. Isosymmetric structural phase transitions: Phenomenology and examples. Acta. 

Cryst. B 51, 753-757 (1995). 
[36] Bleochl, P. E. Projector augmented-wave method. Phys. Rev. B 50, 17953-17979 (1994). 
[37] Monkhorst, H. L. &; Pack, J. D. Special points for Brillouin-zone integrations. Phys. Rev. B 

13, 5188-5192 (1976). 

[38] Sreedhar, K., Honig, J. M., Darwin, M., McElfresh, M., Shand, P. M., Xu, J., Crooker, B. C. 
& Spalek, J. Electronic properties of the metallic perovskite LaNiOa: Correlated behavior of 
3d electrons. Phys. Rev. B 46, 6382-6386 (1992). 



15 



Supplemental Materials 

Experimental results 

The lattice parameters of the LaNiC>3 (LNO) films are determined from x-ray diffraction 
carried out at room temperature. Scans along the (0 L), shown in Fig. [5^,-d, are used to 
determine the c-axis parameters of LNO on LaAlOs (LAO) and SrTi03 (STO). The in-plane 
parameters are determined from HK reciprocal space maps taken about the (2 2 2) film peaks 
(Fig. [5^,f). In both LNO on STO and LAO, the films are found to be coherently strained to 
the substrates, exhibiting the same a,b parameters as the substrates. 

Half-order peaks, arising from the octahedral rotations, are measured to determine the 
oxygen positions. Figure [6] shows scans through various half-order peaks along the H and 
L directions, confirming that the peak widths are independent on their location on the 
reciprocal lattices. 

There are four possible domains that give rise to different structure factors for each 
peak. The domain volume fractions are obtained by comparing the intensities of a set of 
symmetrically equivalent peaks with fixed L. For LNO on STO, the equal intensities of 
the (1/2 3/2 1/2) family of peaks, shown in Fig. [7^,, indicate the four domains have equal 
volume fractions. Similar scans for the LNO/LAO film (Fig. [7^) indicate unequal domain 
populations, with only two domains present in a ratio of 2:1. The ratio of twin domains in 
the LAO substrate is also 2:1 (Fig. [8]) suggesting that the rotational domains present in the 
LNO film are imprinted from the LAO substrate. 

With the domain volume fractions determined, scans were carried out through multiple 
half-order Bragg peaks to determine the oxygen positions. The bond angles and lengths are 
obtained and are summarized in Table SI. 
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Figure 5: X-ray diffraction used to determine the pseudocubic lattice parameters of the LNO 
films. a,b Extended (0 L) scan and data near the (0 2) peak are used to determine the 
c-axis parameter of the LNO film grown on LAO. c,d Extended (0 L) scan and data near 
the (0 2) peak are used to determine the c-axis parameter of the LNO film grown on STO. 
e, Reciprocal space maps taken near the (2 2 1.946) condition confirm that the film grown 
on LAO is strained, f, Reciprocal space maps taken near the (2 2 2.05) condition confirm 
that the film grown on STO is strained. In all cases, the reciprocal lattice units are set to 
those of the substrate. 
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Figure 6: Scans of various half-order peaks through L and H. The peak widths of all the L 
scans (a) are roughly equal, while the widths of all the H scans (b) are roughly equal. From 
this we conclude that the half-order reciprocal lattice points all have approximately the same 
width along H, K, and L. 




Figure 7: Evidence for domains in LNO films. Half-order diffraction peaks from the LNO 
grown on STO (a) and LAO (b). The equal intensity of the (3/2 1/2 l/2)-family of peaks in 
LNO/STO indicate an equal population of tilt domains. The unequal intensity of the (5/2 
3/2 3/2)-family of peaks in LNO/LAO indicate that some domains are preferred over others. 
The (1/2 1/2 5/2)-family of peaks also exhibits unequal intensities. 
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Figure 8: Twinning in LAO substates. a, An optical micrograph of the LNO/LAO sample. 
Twin boundaries run parallel along one of the a-axis of the LAO substrate, b, Schematic view 
of the twin structure of the LAO substrate and the perovskite unit cells, c, The (0 2) perovs j site 
rocking curve of LAO substrate reveals the presence of two dominant orientations . The 
intensity ratio of the two peaks indicates the ratio of the two twin domains is approximately 
1:2. d, The reflection conditions of the twin domains. The rhombohedral tilting angle is 
exaggerated for clarity. 



LNO/STO (XRD) LNO/LAO (XRD) LNO/STO (DFT) LNO/LAO (DFT) 

din-plane (A) 1.968 ± 0.002 1.916 ± 0.005 1.970 1.921 

dout-of-piane (A) 1.933 ± 0.002 1.949 ± 0.002 1.935 1.950 

0in-planeO 165.8 ± 0.5 164.0 ± 2.0 164.7 161.5 

0out-of-piane(°) 159.9 ± 0.6 175.2 ± 0.6 159.3 173.7 

TABLE I: The bond angles and lengths as determined by x-ray diffraction (XRD) and density 
functional theory (DFT). 
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Geometric constraints on the octahedral rotation system 

In Fig. [9^, we show the B0 6 octahedral unit as the building block for the perovskite 
structure. We consider in Fig. [9j3, an ideal cubic structure with an axial ratio c/a denned to 
be unity. In this case all of the B-0 bonds and O-B-0 bonds are equal yielding the ideal 
six-coordinate octahedron. 

In Fig. [9)3 we schematically show the effect of applying bi-axial compressive strain to a 
3x3x3 supercell of the primitive (and cubic) 5-atom perovskite cell. We enforce that the 
B-0 bond lengths and O-B-0 angles in each individual octahedron remain identical. [This 
assumption is good for small strains, and agrees with Pauling's rule on parsimony] In order 
to accommodate the change in the in-plane a lattice constant, the octahedra rotate about 
the c-axis by an amount 7. In the out-of-plane direction, if the octahedra do not change 
their rotational pattern, then by the requirement that the out-of-plane lattice constant c 
increases (from elastic theory), there results a mismatch 5c between the repeat unit cell 
corresponding to c/a=l (dashed line), and that with c/a > 1 (yellow line). This mismatch 
that occurs when the octahedra along c maintain their original rotation (in this case zero) 
and bond lengths, requires that the magnitude of the octahedral rotation pattern about an 
axis parallel to the plane of the strain changes to make the repeat unit cells coincide with one 
another. The need to change the octahedral bond angles, then leads to the abrupt, first-order 
phase transition. Interestingly the origin for this isosymmetric constraint can be completely 
attributed to the corner connectivity of the octahedral units and their antiferrodistortive 
nature that is compatible with a three-fold axis but incompatible with the bi-axial constraint 
imposed on the film by the substrate. 
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Figure 9: Geometric illustration for the origin of the isosymmetric phase transition at the 
axial ratio of c/a—1. a, The BOg octahedral building block in perovskite oxides, b, Ideal 
reference structure for a perovskite with cubic symmetry. The ideal, cornershared octahedral 
units, are shown to readily tile the unit cell in the x and z-directions. c, The effect of 
bi-axial compressive strain causes the a-axis parameter to shrink, and therefore increases the 
c parameter due to elasticity. If we fix the B-0 bond lengths about their ideal length as in 
(b), then the octahedral must rotate by an amount 7 along the ^-direction. The geometric 
mismatch Sc between the ideal and elongated unit cells then drives the octahedral to rotate 
about an in-plane axis, in order to preserve the original B-0 bond lengths and O-B-0 internal 
angles. 
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Figure 10: Calculated electronic structure. The calculated band structure of LaNi03 on 
LaAlOs (a), bulk LaNiC>3 (b), and LaNiC>3 on SrTi03 (c) along high symmetry lines in 
the Brillouin given as r(0,0,0)^M(l/2,l/2,0)-^X(l/2,0,0)^r(0,0,0). For each strain state, 
LaNiC*3 remains metallic. The role of strain on the electronic structure is to partially split 
the orbital degeneracy between the majority spin d x 2— y 2 and d 3z 2— r 2 orbitals (blue points) 
as indicated at the T-point. The Fermi level is denoted by the solid red line at eV. 

Electronic band structure calculations 



We next plot the electronic band structures for the two strained films on STO and LAO 
using the experimental nickelate structures along the high symmetry lines in a non-primitive 
(2x2x1) unit cell of the monoclinic structure. This supercell doubles the number of bands, 
however, it makes the identification of the orbital contributions to the band structure more 



transparent. In Fig. 10 1 we show the calculated band structure for the experimentally 
determined nickelate film on the LAO substrate and compare it to the bulk experimental 
structure in Fig. [TOj j and to that on STO substrate in Fig. [I~0p . For each strain state the 
band structure is found to be fully metallic, with strain partially lifting the majority spin e g 
orbital degeneracy at the T-point around the Fermi level (set to zero here). The metallic 



nature of the films is evident in the resistivity behavior, shown in Fig. 11 
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Figure 11: Resistivity of the LNO films on LAO and STO. Both films show metallic behavior 
with a similar temperature dependence. 
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